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Abstract 

We compare theoretical techniques utilized for description of optical response in molecular junc¬ 
tions, and their application to simulate Raman spectroscopy in such systems. Strong and weak 
sides of the Hilbert vs. Liouville space, as well as quasiparticles vs. many-body states, formula¬ 
tions are discussed. Common origins of the methodologies and different approximations utilized in 
different formulations are identified. 
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I. INTRODUCTION 


The interaction of light with molecules is an important field of research due to its ability 
to provide information on molecular structure and dynamics, and to serve as a control tool 
for intra-molecular processes. Theory of molecular optical spectroscopy has been developed 
and widely utilized in studies of optical response of molecules either in the gas phase or 
chemisorbed on surfaces |1]. 

Recent progress in nano fabrication made it possible to perform optical experiments on 
molecular conduction junctions. In particular, current induced fluorescence |2j and Raman 
measurements [3H5] were reported in the literature. Theoretical description of optical re¬ 
sponse of such nonequilibrium open molecular systems is challenging due to necessity to 
account for optical excitations in the system in the presence of electron flux through the 
junction. For example, correlation between Stokes signal and conductance were measured|3J 
and an attempt of theoretical explanation op the effect was proposed[6j [7]. With scatter¬ 
ing theory being inapplicable to the many-body electronic problem of quantum transport 
in junctions, the very definition of optical scattering processes in such systems is a chal¬ 
lenge. Also, as discussed below, formal application of the molecular spectroscopy theory to 
junctions may be problematic. 

Here we compare theoretical approaches utilized in the literature for description of op¬ 
tical spectroscopy of current carrying molecular junction. After introducing model of the 
junction subjected to external radiation field in Section [TT] we consider theoretical founda¬ 


tions of molecular spectroscopy in Section |III| We then discuss Hilbert vs. Liouvillc space 
formulations in Section [TV] and quasiparticles vs. many-body states versions of the two for¬ 
mulations in Section [VJ In Section |VI| we specifically focus on different contributions to the 
total optical signal, and identify those relevant for the Raman scattering. Our conclusions 
are presented in Section VTTJ 


II. MODEL 

We consider a junction formed by a molecule M coupled to two contacts L and R (each 
at its own thermal equilibrium) subjected to an external radiation field. Hamiltonian of the 
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model is 


H =H 0 + V 

(1) 

H 0 =H S + V (h b + V sb ) + H rad 

(2) 

B=L,R ' ' 


v = ( U M, a O j M a a + H.c.) 

(3) 


M,a 


where H 0 describes molecular junction and radiation field, and V introduces coupling be¬ 
tween them. Here H $, Hb , and H ra d are Hamiltonians of the molecule (system), contacts 
(baths), and radiation field, respectively. Hs describes electronic and vibrational structure 
of the molecule as well as arbitrary intra-molecular interactions. The contacts and field are 
assumed to be reservoirs of free electrons and photons, respectively 


Hb=^2, £ kc\c k 

(4) 

k&B 


H ra d ^ ^ ^ a^ofla. 

(5) 


a 


Vsb couples molecule with contacts and usually is assumed to be quadratic in the quasipar¬ 
ticle representation 

Vsb ^ f V km c ^d m T V mk djy, c k j (6) 

meS,keB ^ ' 

In Eqs. |l])-([6j d) m (d m ) and c\. (c k ) are the creation (annihilation) operators of electron in 
level m of the molecule and state k of the contact, respectively. a) a ( a a ) creates (destroys) 
photon in mode a of the radiation field, 0' M ( Om ) is operator creating (destroying) molecular 
optical excitation M, and Um,o is coupling between the excitation M and mode a of the 
radiation field. Here and below h — e — 1. 


III. OPTICAL SIGNAL AS PHOTON FLUX 

Similar to considerations of electron transport in junctions, theories of optical spec¬ 
troscopy are mostly focused on simulating fluxes. More advanced treatments involve also 
higher cumulants of the counting statistics [8]. Diagrams for photon (boson) fluxes are pre¬ 
sented in Fig. [lj As usual flux is defined as rate of population change in the continuum of 
modes {a} of the radiation field, and corresponding expression for the photon flux from the 
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FIG. 1: Photon flux diagrams in (a) quasiparticle and (b) many-body states representations. 
Solid and wavy lines represent electron and photon propagators, respectively. Non-directed lines 
indicate both possible directions, t is the time of the flux, t' is the Keldysh contour integration 
variable. Summation over all indices and integration over contour variables is assumed for every 
connection in the diagrams (except the fixed time of the flux t). 


system into modes of the field at time t is PHD] 


d 


J (t) = ^j t ^ a (t)a a (t)) 


=2Re / dt'Tx r>u 
J — OO 




(7) 


Here Tr[. ..] is trace over molecular optical excitations M (depending on the formulation 
these are either of quasiparticle type or represent transitions between many-body states 
of the molecule), n < ^ > ^ and G < ( > ' ) are lesser (greater) projections of the self-energy due 
to coupling to the radiation held and correlation (Green) function of molecular optical 
excitations, respectively m 


nMi,M 2 ( r ? Tji(*T , t) Ucx^M2 

a 

Gm u m 2 {t, t’) = - i(T c d Ml ( t) d\ h (T')) 


( 8 ) 

(9) 


where r and r' are the Keldysh contour variables corresponding to real times t and t', T c is 
the contour ordering operator, and 


F a (r , ,T) = -i(T c d a {r') a^(r)) 


( 10 ) 
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is the Green function of free photon evolution. Its lesser and greater projections are 


( 11 ) 

( 12 ) 


F<(t',t) = -iN a e* ,a V- t) 

) = -i[l + iY Q ] e -^G-d 


where N a is average population of the mode a. 

We note that expression (J7]) is exact. It consists of two contributions: in-scattering 
(photon absorbed by the system; first term in the second row) and out-scattering (pho¬ 
ton emitted by the system; second term in the second row). The former characterizes e.g. 
absorption spectrum, while the latter can give information on fluorescence. Standard treat¬ 
ment proceeds by evaluating the molecular correlation functions Q < ^ > - > in the presence of 
intra-molecular interactions, coupling to environment (e.g. molecular coupling to contacts in 
junctions), and radiation held. Depending on the formulation some or all of these processes 
are taken into account in the molecular Green function through corresponding self-energies 
E. The latter usually can be evaluated only approximately. Corresponding expressions 
should be derived from the Luttinger-Ward functional na. so that the resulting approxi¬ 
mation fulfils conservation laws pa mj. The molecular Green function can be evaluated 
utilizing the Dyson equation 


G(t,t') = Gn{r, t') + dri dT 2 Go(T,r 1 )'E(T 1 ,T2)G(T2,r / ) 


(13) 


where Go is the Green function in the absence of interactions. In most cases Eq. (13) has to 
be solved self-consistently due to dependence of the self-energy E on the molecular Green 
function. Substituting the converged result into ([7]) yields information on the incoming, 
outgoing, or total optical flux. 

In practice treatment of molecular correlation function is often done in a simplified manner 
directly employing perturbation theory (at least in coupling to the radiation held). While 
this way conserving character of the approximation cannot be guaranteed (e.g. lowest order 
in electron-photon interaction - the Born approximation - is known to be non-conserving 
and in some cases may be not fully adequate B9). the simplification avoids necessity of 
self-consistent treatment, and allows direct classification of multi-photon processes in the 
system. For example, expanding the Green function Q, Eq. (§■ in pertubation series up to 
the second order in interaction with the radiation held V, Eq. (p|, yields 


G(t,t') = Go(t,t') - dr i / dr 2 n(ri,r 2 ) (T c 6(r)d t (ri)6(r 2 )6 t (r / ))c 


(14) 
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where subscript 0 indicates that corresponding averages are evaluated in the absence of the 
field (intra-molecular interactions and coupling to contacts are still present). Substituting 
first term in the right into ([T]) yields lowest (second) order contributions into absorption 
(in-scattering) and emission (out-scattering) spectrum. Substitution of the second term in 
the right provides information on all fourth order optical processes in the system. 

Below we focus on the out-scattering photon flux, second term in the bottom row of 
Eq. (JtJ) , and discuss fourth order optical processes, second term in the right of Eq. (14), 
in the emission spectrum within the Hilbert and Lioville space formulations. Following the 
tradition we pick one of the modes (or set of modes at particular frequency) in the sum (J8]) , 
substitute it in place of the self-energy in second term in the bottom row of Eq. ([?]), and 
assume the mode (which we call final, /) to be empty. Similarly, in Eq. (14) we pick one of 
the modes (which we call initial, i ) from the sum and use it in place of the self-energy in 
the expression. This leads to 

Jir+f(t) =2 Re [ 


dt' 


dn j dr 2 ^ UM'jU fM UM 1 ,iUi,M 2 Fj(t' -t)F i (r 1 ,T 2 ) 


(TcdUf) 6u(t) dL M Om 2 (t 2 )) 


(15) 


We note in passing that separation of the self-energies n into modes in principle can be 
avoided (see e.g. Refs. USED- 


IV. HILBERT VS. LIOUVILLE SPACE FORMULATION 

Difference between Hilbert vs. Liouvillc space formulations stems from the way projec¬ 
tions of contour variables are performed. The Hilbert space formulation utilizes the Keldysh 
contour, and thus transition from contour variables to real times (i.e. ordering positions of 
the times on the contour) requires to work with several types of projected correlation func¬ 
tions (depending on positions of the times on the contour). The Liouville space formulation 
works with the real time axis, ordering the times relative to each other on the axis, and 
thus necessitates introduction of Liouvillc superoperators, whose role - distinguish between 
action on different branches of the contour - is exactly the same as that of different corre¬ 
lation functions projections in the Hilbert space formulation. Naturally, both formulations 
yield the same result [T81423] . 
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FIG. 2: Contour projections for the fourth order interaction with the modes i (red, times t\ and 
t 2 ) and / (blue, times t and t') of the optical held. Time increases from left to right in these 
diagrams. 


Utilizing the Hilbert space procedure in Eq. (15) is equivalent to consideration of all 
possible placements of times t\ and t 2 (real times corresponding to contour variables T\ and 
T 2 ) between times t and t' on the Keldysh contour. Note that in Eq. (15) the latter are placed 
in such a way that t' follows t on the contour (this results form the greater character of the 
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photon Green function F^(t',t)). These placements lead to 12 possible diagrams (i.e. time 
orderings on the contour) shown in Fig. [2j One can also get 12 additional diagrams, which 
are just complex conjugate versions of those presented in Fig. [2] (they come from the Re in 
Eq. (15)). Causal character of interactions implies that time t (the time of the signal) is the 
latest time, usually placed at the point where the two branches of the contour meet each 
other, t' then belongs to the anti-ordering branch, and t\ and 0 may populate any of the 
three regions (preceding t, between t and t! and after t' on the contour). Also times t and t' as 
well f i and t 2 are connected by lines representing free photon propagation described by Green 
functions F. Character of the Green functions (greater for the pair t and t/, Fj?(t',t), and 
either greater or lesser for the pair t\ and t 2 , F^ < \ti,t 2 )) depends on ordering of the times 
on the contour and represents photon emitted (for the greater projection, out-scattering) or 
absorbed (for the lesser projection, in-scattering) by the system. 


The Liouville space procedure applied to Eq. (15) deals with the same problem of ordering 
variables T\ and r 2 between the two times t and t ’. However this time the ordering is 
performed along the real time axis (i.e. not only relative position of times on the contour 
but also relative position on the real time axis is tracked), thus number of diagrams (different 
orderings) is bigger here. Information on the branch of the contour is provided by the two- 
side Feynman diagrams (the Keldysh contour diagrams with additional restriction on the 
relative positions of times with respect to the real time axis). It is customary to indicate 
each photon process by separate arrow in these diagrams, rather than consider contractions 
representing free photon propagation. The agreement is that arrow pointing to the left 
corresponds to creation operator of the photon in quantum mechanical description of the 
field (or factor e lut for classical treatment of the field), while arrow pointing to the right 
represents operator of annihilation of the photon (or factor e~ wt ) [Tj. The procedure being 
applied to expression (15) leads to 24 double-sided Feynman diagrams presented in Fig. [3j 
These are the diagrams presented in Figs. 4 and 5 of Ref. [2TJ Similar to the Hilbert space 
formulation 24 more diagrams are complex conjugate variants. 

Comparison of the two sets of diagrams is straightforward, when one keeps in mind an 
additional restriction on time ordering within the Liouville space formulation (ordering not 
only on the contour but also with respect to real time axis). So, one Hilbert space type 
diagram (see Fig. [2]) may encompass several diagrams of the Liouville type (see Fig. [3j) . 
Diagram-to-diagram correspondence is shown in Fig. |4j 






FIG. 3: Contour projections (double-sided Feynman diagrams) for the fourth order interaction 
with the modes i (red, times t\ and U) and / (blue, times t and t') of the optical field. Time 
increases from bottom to top. 


V. QUASIPARTICLES VS. MANY-BODY STATES APPROACH 


In realistic simulations in junctions in addition to molecular interaction with radiation 
field other intra-molecular interactions (e.g. electron-electron or electron-phonon) and cou¬ 
plings to environment (e.g. coupling to contacts) should be accounted for. 

Quasiparticles approaches (in either Hilbert or Liouville space) are usually efficient in 
treating molecule-contacts couplings. Indeed, since the latter are represented as quadratic 
in terms of elementary excitations, Eq. they do not alter non-interacting character of the 
Hamiltonian, and thus Wick’s theorem is available to exactly evaluate multi-time correlation 
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FIG. 4: Comparison between contour projections presented within the Hilbert (see Fig. [ 2 ]) and 
Liouville (see Fig. [ 3 ]) space approaches. Time increases from bottom to top. 

functions in presence of the couplings [25] . However, quasiparticles are less convenient 
when dealing with intra-molccular interactions. That is, if intra-molecular interactions are 
weak compared to either molecule-contacts couplings or separation between electron energy 
from molecular resonances, standard diagrammatic technique (on the contour perturbation 
theory) can be invoked to account approximately for the interactions. 

For strong interactions (e.g. Coulomb blockade or polaron formation) usually encountered 
in the resonant tunneling regime, a nonequilibirum atomic limit (i.e. consideration utilizing 
many-body states of the isolated molecule as a basis) is preferable [2BJ. Note that historically 
spectroscopy of isolated molecules also was formulated mostly in the language of many- 
body states (molecular or dressed states of molecule and the field) p]. Many-body states 
based formulations, while accounting exactly for intra-molecular interactions, are capable of 
treating molecule-contacts couplings only approximately. We note that such approximate 
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schemes exist for both Hilbert [27l 129] and Liouvillc [301 [32] space formulations, however 
these approximations are not well controlled. 

A promising tool is the Hilbert space type pseudoparticle NEGF (PP-NEGF) tech¬ 
nique [331 i36] , which while being a many-body states formulation allows to take into account 
coupling to the contacts within well controllable diagrammatic technique, which is based on 
extension of the standard methods of quantum field theory into extended version of the 
Hilbert space. The PP-NEGF has several important advantages: 1. The method is concep¬ 
tually simple; 2. Its practical implementations rely on a set of controlled approximations 
(standard diagrammatic perturbation theory techniques can be applied); 3. Already in its 
simplest implementation, the non-crossing approximation (NCA), the pseudoparticle NEGF 
goes beyond standard QME approaches by accounting for both non-Markovian effects and 
hybridization of molecular states; 4. The method is capable of treating junction problems 
in the language of many-body states of the isolated molecule, exactly accounting for all the 
on-the-molecule interactions. Recently a first application of the methodology to problems 
of optical spectroscopy of junctions was proposed in Ref. [37, 


VI. RAMAN SPECTROSCOPY 


Eq. (15) represents all contributions to the fourth order optical process. Note that in¬ 


teraction with mode / only includes outgoing photon, while interactions with mode i have 
both incoming and outgoing photon contributions. Correspondingly, the projections can be 
separated into those containing two outgoing photons (diagrams B in Fig. [2j) and those which 
contain incoming photon of mode i and outgoing photon of mode / (diagrams A in Fig. [2j) . 
The former contribute to fluorescence (both coherent fourth order and sequential second 
order processes), while the latter yield Raman (see below) and absorption with emission 
(sequential second order) processes. 

Both diagrams A and B of Fig. [2] contain virtual photon processes. For example, interac¬ 
tion with mode i in projections (a4) (A type diagram) and (al) (B type diagram) of Fig. [3] 
are virtual photon processes. These processes yield renormalization of molecular correlation 
functions due to presence of the radiation field. Accounting for such processes within the 
procedure described above should be done with care. Indeed, reasonable estimate of the in¬ 
teraction with the field is U ~ 1CT 3 — 10 -2 eV 1331. For a molecule chemisorbed on metallic 
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surface electronic escape rate (characteristic strength of the molecule-contact interaction) is 
T ~ 0.01 — 0.1 eV [39], and in this case renormalizations due to coupling to contacts are 
much more pronounced. For a molecule detached from contacts renormalizations due to the 
radiation held become important, however proper way to account for those is by solving 


Eq. (13). Accouting for the renormalization by keeping only second term in the right of 


Eq. (14) (besides being a non-conserving approximation) results in a Green function diver¬ 
gent at molecular resonances. Resummation of the perturbation series and introduction of 
the retarded projection of the self-energy is the proper way to account for the renormaliza¬ 
tion, when coupling to the radiation held is weak relative to intra-molecular interactions. 
If coupling to the radiation held is the dominating interaction, perturbation series treat¬ 
ment becomes invalid. We stress that in experiment one deals with the total optical signal, 
and a proper theoretical treatment should utilize Eq. ([7]) with conserving approximations 
employed to evaluate molecular Green functions. 

We now focus on theoretical description of a particular optical process - the Raman 
scattering. Raman scattering from mode i to mode / is a coherent process of (at least) 
fourth order in coupling to radiation held with two orders in coupling to incoming photon 
of mode i and two orders in coupling to outgoing photon of mode /. The process should 
also satisfy a proper energy conservation: at steady-state difference between incoming and 
outgoing photon frequencies should be equal to multiples of the frequency of molecular 
vibration (may be plus electronic energy differences within the same electronic level or 
many-body state broadened by coupling to environment) 


jRaman ~ ^ + ^ + 


(16) 


ffere n is integer number (n — 0 yields Rayleigh scattering) and A E is electronic energy 
change within the same broadened level or many-body state. 

Identification of diagramms related to Raman scattering requires some care. An easy 
way to identify relevant contributions relies on the Hilbert space formulation of section |IV| 
with utilization of the Langreth projection rules m First, only projections A of Fig. [2] 
are those corresponding to a process with one incoming and one outgoing photon. Second, 
among those only parts of projections A2, A3, A4, and A5 yield proper energy conservation, 


Eq. (16). For example, let employ the Langreth rules to projection A4. This yields 4 terms 
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in expression (15) 


r+oo 

' — oo 
/»+oo 

' — oo 
/*+oo 


' —oo 
/»+oo 


+ 


/ 0 /*0 

G?(ti ~t) d(t 2 - h)... 

oo J —oo 

/ 0 /*0 

d(ti — f) / d(t 2 ~ t')... 

oo J —oo 

/ 0 pO 

d(ti — t! ) / d(t 2 — ti)... 

-OO J — OO 

/ 0 /»0 

d(ti — / d(t 2 — t')... 

oo J —oo 


(17a) 

(17b) 

(17c) 

(17d) 


Only (17b) yields proper energy conservation. Similarly, application of the Langreth rules 


to projection A3 yields one, and to projections A2 and A5 yields two integrals for each with 
only one term in each case satisfying restrictions of the Raman process. 

In Refs. HD], HO] we presented a theory of Raman scattering in molecular junctions formu¬ 
lated within the NEGF (the Hilbert space quasiparticle formulation). Projections presented 
there were different from those given in Fig. [2] in two aspects: a. the projections there are 
given after the Langreth rules have been applied and b. some of the projections (see e.g. 
Fig. 8b of Ref. HOl) represent ‘the hole view’ of scattering process. We note that in general 
(when explicit time-dependent processes or a magnetic held are present) the latter flexibil¬ 
ity does not exist. However for steady-state situation the two pictures can be shown to be 


equivalent. For example, the contribution (17b) to the Raman signal is (note, it is easy to 
see that the contribution is real following the derivation leading form Eq. (28) to Eq. (56) in 

Ref. rimi 


I UiU f \ 2 Ni 


^+00 rO rO 

d(t'-t) / d(t!-t) / d(t 2 - t') 

-OO J—OO J—OO 

X (6t(P)6(t 2 )6t (tl) o (t) ) 

r»+00 


= | Ui U f \ 2 Ni / d(t’ - t ) / d(t i - t) / d(t 2 - t') e-fc'R*-*') e ^-t2) 


X (&(t)d(t 1 )d\t 2 )d(t’)) 


where N t is the population of the mode i. Then utilizing the time-reversal symmetry of the 
molecular correlation function. |IT| 

( 0 \t) 0 (h) o\t 2 ) 6(f)) = << 6\-t’) d(-t 2 ) o\-ti) 0(-t)) (is) 

and inverting signs of the time variables one gets Eq. (29) of Ref. (10], which was derived 
from ‘the hole type’ projection shown in Fig. 8b there. We note that while the results are 
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equivalent, graphical forms of the projections differ: Fig. 8b of Ref. [TO] has the two incoming 
field modes on different branches, while the diagram A4 - on the same branch. 


VII. CONCLUSION 

We discussed theoretical formulations utilized in studies of optical processes in current 
carrying molecular junctions. We start by identifying total optical signal as photon flux from 
the system to registering device, for which exact formulation in terms of molecular correlation 
(Green) functions is available. Standard quantum field theory methods allow to formulate 
well controlled conserving approximation to evaluate the Green functions of the molecule. At 
the same time tradition coming from the theory of optical spectroscopy for isolated molecules 
utilizes perturbation theory. The latter, while not being a conserving approximation, allows 
to separate different photon scattering events in the total optical signal. We discussed the 
Hilbert and Liouvillc space formulations, indicating similarities and differences in graphical 
representations of photon scattering process in the two formulations. We also pointed out 
strong and weak sides of quasiparticles (second quantization) vs. many-body states based 
(nonequilibrium atomic limit) approaches in both formulations, and outlined the PP-NEGF 
as a promising theoretical approach to study optical spectroscopy in junctions. Finally, we 
focused on theoretical formulation for Raman scattering in molecular junctions, and clarified 
questions raised in the literature m with respect to identification of relevant projections and 
obtaining corresponding contributions to the total optical signal. Thus, the paper bridges 
different formulations of optical spectroscopy in open nonequilibrium systems. 
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